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Abstract 

We define generalized Husimi distributions using generalized coherent states, and show that 
their moments are good measures of complexity of many-body quantum states. Our construction 
of the coherent states is based on the single-particle transformation group of the system. Then 
the coherent states are independent-particle states, and, at the same time, the most localized 
states in the Husimi representation. Therefore delocalization of the Husimi distribution, which 
can be measured by the moments, is a sign of many-body correlation (entanglement). Since the 
delocalization of the Husimi distribution is also related to chaoticity of the dynamics, it suggests a 
relation between entanglement and chaos. Our definition of the Husimi distribution can be applied 
not only to the systems of distinguishable particles, but also to those of identical particles, i.e., 
fermions and bosons. We derive an algebraic formula to evaluate the moments of the Husimi 
distribution. 

PACS numbers: 03.67.Mn, 05.45.Mt 
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I. INTRODUCTION 



Although the manifestation of chaos in quantum systems has been investigated exten- 
sively in the past few decades, definition of quantum chaos is still unclear. Especially the 
studies of many-body systems are far behind those of one-body systems. 

In 0, we proposed a measure of complexity (chaoticity) of quantum states in one-body 
systems by using the Husimi distribution 2f . The purpose of this paper is to generalize this 
idea to many-body problems. We will show that the generalized Husimi distribution is also 
related to quantum many-body correlation (entanglement). Therefore it implies a relation 
between entanglement and chaoticity of the dynamics. 

n one-body systems, the Husimi distribution function of a quantum state \(p) is defined 



as 

H M (p,q) = \(p,qM 2 , (1) 

where \p, q) is a coherent state whose average momentum and coordinate are respectively p 
and q. Since chaoticity of classical systems can be characterized by derealization of trajec- 
tories, derealization of the Husimi distribution can be regarded as a quantum manifestation 
of chaos. We used the second moment as a simple measure of the derealization. The inverse 
of the second moment represents the effective volume occupied by the Husimi distribution, 
and has a good correspondence with chaoticity of the classical system. 

To extend this idea to many-body systems is not difficult, at least formally. If we define 
coherent states in a many-body system, the Husimi distribution is defined as the square of 
the absolute value of the coherent state representation. There are several ways to generalize 
the idea of coherent state to many-body systems. We construct the coherent states based 
on the single-particle transformation group, following the group-theoretical construction by 
Perelomov jjj. Then the coherent states are independent-particle states, as we will explain 
below. 

Let us consider a system of m qubits as an example. In this case, the single-particle 



(local unitary) transformation group is S U(2) x • • ■ x SU(2). Coherent states are generated 
by applying this group to the "vacuum" |0) (g> • • • <8> |0). Then we obtain all separable 
(disentangled) states. Hence, 

coherent state •<=>■ disentangled state. (2) 



FIG. 1: The Husimi distribution of qubits is denned on S 2 x S 2 x • • • x S 2 . A disentangled state 
is represented by a point on this manifold, and the Husimi distribution thereof is localized around 
the point. An entangled state is represented by a delocalized distribution. 

The coherent states are parametrized as 

|0» = |01,0l)®---®|0m,0m), (3) 

where \9,4>) is the Bloch sphere representation of a qubit. The Husimi distribution of a 
quantum state \(p) is defined as 

Hm(o,4>) = \(o,4>\<p)\ 2 (4) 

on the "phase space" S 2 x ■ ■ • x S 2 . A disentangled state is represented by a localized wave 

packet in the phase space. (Fig. [TJ) 

The Husimi distribution can also be defined for systems of identical particles, i.e., bosons 

and fermions. Suppose there are m identical particles in N single-particle states. In this 

case, we cannot operate each particle separately. Therefore the single-particle transformation 

™ I — I 

group of this system is not SU(N) x ■ • ■ x SU(N), but SU(N) (or U(N)) [ig|. 
states are defined as 



The coherent 



IO = tf(C)|o>, (5) 

where U(() is an element of U(N) specified by the parameter (. The "vacuum" |0) can be 
written explicitly as 

, nN j \<Pi)\<Pi)---\<Pi) (boson) 

|0) = < , (6) 

I A{\(pi)\<f 2 ) ■■■\ < fm)) (fermion) 

where \(fi) is the i-th single-particle state and A is the anti-symmetrization operator. Then 
the coherent state is 

|c) , |^i(C))l^(0)... 1^(0) (boson) 

^ ^(ki(O)l^(C)) •••km(C))) (fermion) ' 



where |<£>i(0) = ^(Olv 9 *)- I* 1 the bosonic case, the coherent states are separable, and it is 
easy to see that all separable bosonic states can be written in this form. Therefore 

coherent state separable state. (8) 

In the fermionic case, a coherent state is a Slater determinant, and any Slater determinant 
can be written coherent state. Namely, 

coherent state -<=>• Slater determinant. (9) 

In any case, the coherent states are the least correlated ones. We summarize the features of 
these systems in Table. H] 

The Husimi distribution is defined by using the coherent state as 

n M (() = \mf . (io) 

Since the most localized states in the Husimi representation are the coherent states, dereal- 
ization of the Husimi distribution implies correlation among the particles. The derealization 
can be measured by the Renyi-Wehrl entropy |4( 

where Afw is the moment with an index q > 0: 

M$ = J dn(C) {H M (0Y ■ (12) 

Here, dfi(Q is the Haar measure of the group. Hence is an invariant of the group 

transformation. Note that reproduces the normal entropy in the limit q — > 1: 



Jd^)H W) {Q\nH W) {Q. (13) 



The Renyi-Wehrl entropy represents the effective volume occupied by the Husimi distri- 
bution. For instance, if the Husimi distribution takes the same value over a region with 
volume V and takes zero value outside of it 20] , = \nV. Therefore it is natural to 
expect that takes the minimum value for the coherent states. This conjecture (gener- 



alized Lie 
and Lieb 



)-Wehrl conjecture) was formed in 4], following the questions raised by Wherl 

7[. According to the conjecture, 



For integer indices q > 2 it was proved in 

, n s r „s separable (qubits, bosons) 

S$ = S { ± <=► \<p) : { W ' (14) 

I Slater determinant (fermions) 



Thus S^ q ' shows whether the particles are correlated or not. If is large, the Husimi 
distribution is delocalized, and we need a lot of wave packets (independent-particle states) 
to represent the state. Threfore the Renyi-Wehrl entropy can be regarded as a measure of 
complexity. 

Note that the complexity we measure here is completely independent of the complexity 
of single-particle states which was the subject of [jj. For instance, we can construct a Slater 
determinant from highly chaotic single-particle states, which is still the most localized state 
in the many-body Husimi distribution. 

When a Hamiltonian of the system is specified, the time evolution operator of the system 
can be represented by the coherent state path integral. The stationary phase condition of 
this integral leads to an equation of motion in the phase space, which defines the "classical" 
dynamics of the system. (See, for example, JfJ.) In the case of qubits, for example, the 
equation of motion is that of classical spins. When the "classical" motion is chaotic, the 
Husimi distribution will probably spread all over the phase space. This means that chaotic 
dynamics leads to highly complex states. 

Next we consider how to calculate the moments (|12j) . Although the definition of the 
moments (|12J) using the Haar measure is an elegant way to obtain invariants of the group, 
this integral representation is almost useless in numerical calculations because the dimension 
of the phase space is huge when there are many degrees of freedom. Since the coherent states 
form an overcomplete set, the coherent state representation is highly redundant. We should 
avoid such redundancy in real calculations, as we did in Q| for one-body problems. 

A complex quantum state is usually represented by expansion coefficients in a basis. 
Therefore it is desirable to represent the moments (fT2|) directly by them. Fortunately, it 
is possible if the index q is a positive integer. We will derive an algebraic formula which 
represents the moment directly by the expansion coefficients in section |H] This explicit 
formula is not only practically useful, but also theoretically insightful. 

This paper is organized as follows. In section |H] we derive the explicit formula for the 
moments. In this formula, the moments are represented by the expansion coefficients of the 
state and group-theoretical factors which are independent of the state. The main concern 
in the following sections is to determine the group-theoretical factors. After explaining the 
general formalism, we consider bosonic systems with two single-particle states in detail in 
section ITTT1 These systems have U(2) ~ {7(1) x SU(2) as the single-particle transformation 



5 





one-body system 


bosons 


fermions 


distinguishable particles 




(n dim.) 


B[N, to] 


F[N, to] 


D[N U . 


..,N m ] 


group 


HW group 


U(N) 


U(N) 


SU(Nx) x •• 


■ x SU(N m ) 


phase space 




C pN-l 




CP* 1 " 1 x •• 


• x CP Nm - x 


coherent state 


Gaussian 


separable state 


Slater Det. 


separable state 




wave packet 


\<p)\<p)...\tp) 


<A(\<Pl)\<P2)---\<Pm)) 


lvi)|ya> 


■ ■ ■ \<Pm) 



TABLE I: Coherent states for one-body systems, bosons, fermions and distinguishable parti- 
cles. B[N,m] (F[N,m\) means a system of to bosons (fermions) in N single-particle states. 
D[N%, . . . ,N m ] is a system of to distinguishable particles in which i-th particles takes iVj states. 
CP n is the n-dimensional complex projective space and Cr mi jv(C) ~ U(N)/ (U(N — to) x U{m)) 
is the Grassmann manifold. Note that CP 1 ~ S 2 . HW group in the second row means the 
Heisenberg-Weyl group, which is generated by p and q. 

group, and the corresponding phase space is the two-demensionasl sphere S 2 . Therefore we 
can visualize the Husimi distribution in this case, which will help us understand the idea 
of the Husimi distribution in many-body systems. Based on this analysis, we treat bosonic 
systems generally in section IIVI In section IS we treat fermionic systems. In section IVI1 we 
consider distinguishable particles including qubits. We investigate the second moment in 2- 
and 3-qubit cases in detail, and show the relation between the second moment and other 
known measures of entanglement. The final section is devoted to concluding remarks. 



II. EXPLICIT FORMULA FOR THE MOMENTS WITH INTEGER INDICES 

In this section, we show how to obtain an explicit formula for the moment for a 

positive integer q. Our method is based on the expression of the moments ([21)1 using the 
tensor products of the state, which was used also in [t| to prove the generalized Lieb-Wehrl 
conjecture. We briefly review the derivation of (|2ip. and show how the formula is obtained 
from it. 

To begin with, we briefly state the definition of the generalized coherent states of a 
compact semisimple Lie group G [3]. The Lie algebra of G can be written in Cartan basis 
{E a ,Hj}. Let D\ be an irreducible representation space of G with the lowest weight —A. 
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Coherent states in D\ are obtained by the action of the group G on the lowest weight state 
| — A), which can be written explicitly as 

|C,A)=W(Oexp(C a £ a )|-A). (15) 

Here, Af(C) is the normalization constant and a runs over all positive roots. 
Next we consider the moment with a positive integer index q 

M$=c x j dfi(C) |(C,A|^)| 9 . (16) 

Here, dfi(0 is the Haar measure of G and c\ = dimD^ is the normalization constant. The 
moment can be rewritten as 

M® =c x j dMC) W m ) (IC, A)^) ((C, AD (<p)®«) . (17) 

The important point to note is that the tensor product of the coherent state IC)® 9 is again 
a coherent in another irreducible representation D q \ C whose lowest weight is — qX. It 
can be shown explicitly as 

IC, Ar = Af(C) 9 exp (( a E^) | - A}® 9 = |C, gA). (18) 

Here, is the representation of E a in -D® 9 , whose explicit form is 

E^ = E a ® J ® • ■ • ® J + / ® £ a ® • • • ® J H h I ® • • • ® I ® £«. (19) 

It is easy to verify that | — A)® 9 is the lowest weight state of D q \. 
Since the "resolution of unity" 

/ = c gA y dMC)IC,<?A)(C,gA| (20) 

holds in D q \, (JT7j) can be written as 

<'> = aErz£ ^ Pd <> ^ • (21) 

where Pd qX is the projection operator to D q \. For our purpose, it is convenient to introduce 
normalized moments 

~ () ^dnn^ () = 
l v > dimD A M v VqXh K J 
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where the symbol () denotes the expectation value for ly?)® 9 . Then = 1 if \<p) is a 

coherent state, and M$ < 1 if it is not. 
If the state is expressed as 

|¥>> = X>|i) (23) 

i 

with an orthonormal basis {\i)} in D\, the normalized g-th moment is 

M,S = El^| 2 ' ( 24 ) 



b j = E 0'? A Ni'---'^) c n ( 25 ) 

il,...,iq 

where {\j,q\)} is an orthonormal basis of D q \, and 

(j, q\\h,..., i q ) = (j, qX\ (|zi) ® |z 2 ) <g> • • • ® |z 9 )) . (26) 

Thus we have obtained the explicit formula for the moments. 

Since the terms in the sum (125(1 are invariant under permutations of ix,i2, ■ ■ ■ ,i q , the 
same term appears q\ times if all indices are different. In general, the same indices appears 
in the set {i s }- We represent it as as (i* 1 ,^ 2 , . . . ,zf ! ), which means the index i\ appears k\ 
times. Then this term appears ( fei k q k J times in the sum (J25J) . where 



9 \ _ I kl \ki... kl \ (if fci + fc 2 + • • ■ + fc t = g) 
fei, Afc, • • • , h) ) (else) 



is the multinomial coefficient. Therefore (l25l) can be rewritten as 



(27) 



B i = E ( K k l . . , k ) & 4 2 , • • • , . . . <* (28) 

where 

(j, qX\i1\ ...,#) = (j, qX\ (Izx)^ ® |z 2 )^ 2 ® • • ■ ® . (29) 

The sum in (|28|) is taken over all possible combinations of {i s } and {k s }. 
Let the irreducible decomposition of D® q be 

D T = J2 D ^- ( 3 °) 



V.T 



Here, D VjTv is an irreducible representation with the lowest weight —u, and the additional 
integer suffix r v is the multiplicity label. Then 

y Z{Pu,r,)=J^ q (31) 

holds, where P VyTv is the projection operator to D VjTv and M = (£\ IqI 2 ) 1 ^ 2 is the norm of 
\ip) which is usually normalized to 1. Hence 

M\$) = l ~ E ( p ^>- ( 32 ) 

Note that D q \ is multiplicity-free. In many cases, this formula is more useful and informative 
than calculating (P q \) directly, as we will see later. 



III. SIMPLE EXAMPLE: BOSONS IN TWO SINGLE-PARTICLE STATES 



A. Algebra of generators 

Let us consider a system of m bosons in two single-particle states as a simple example. 
In this case, there are m + 1 many-body states, which are |0,m), |l,m — 1),..., \m,0) in 
the occupation number representation. The many-body states are generated by applying 
particle-hole excitation operators 

Xl = a i a) (1<M,<2) (33) 

to the "vacuum" |0,m). Since the total number operator TV = a\ai + a\a,2 corresponds to 
the irrelevant total phase, we can remove it from the algebra. The rest forms the Lie algebra 
of SU{2): 

J+ = a ± al, (34) 
J_ = a2d\, (35) 

J z = ^(o 2 4 - a i a l)- ( 36 ) 
Actually, it is easy to verify the following commutation relations 

[J Z ,J±] = ±J±, [J + ,J-] = 2J Z . (37) 
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The two single-particle states |0, 1) and |1,0) form the fundamental (spin 1/2) represen- 
tation of SU(2). The m-particle states form an irreducible representation with spin m/2, 
which we denote as -D m / 2 - Hereafter, we mainly use spin quantum numbers j and j z , instead 
of occupation numbers, to specify a quantum state. We put suffixes n and s to distinguish 
the two notations. The relation between the two notations is 



\nx,n 2 )r 



n 1 + n 2 . ni - n 2 x 

3 = o ,3z= o / • ( 38 ) 



B. Coherent state and the Husimi distribution 

The SU(2) coherent states in the spin j representation are defined as 

|j,C)=Ar(C)exp(CJ + )|j,-j) s , (39) 

where A/"(C) is a normalization factor. In our case, the total spin j is determined by the 
number of particles by j = m/2. 

( is considered to be a coordinate system of CP 1 ~ S 2 , which is related to the angular 
variables of the sphere (9, if) by 

C = -e-^tan^. (40) 
The coherent state \j, () is expanded as 

li.O = E /./<).,. (4i) 

t*=-3 

where 

uUO = ( , 2j _. ) (- sin t) J+ " ( cos ?V . (42) 



J + /V v 2; V 2 

The Husimi distribution function of a m-particle state \<p) is defined as 

H W) {Q = \(j=m/2XW- (43) 

For example, let us consider the simplest nontrivial case m = 2 which corresponds to the 
spin 1 representation. The Husimi distributions of the three basis states are 

n mn (e,4) = sin 4 ^, (44) 

W|i,i)„(M) = 2 sin 2 ^ cos 2 ^ = isin 2 #, (45) 

H Mn (B,<l>) = 00s 46 -, (46) 
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10,2 > II, 1> 12,0 > 



FIG. 2: Rough sketch of the Husimi distribution for bosons with N = m = 2. |2, 0)„ and |0, 2) n 
are coherent states, but |l,l) n is not. Therefore |l,l) n has a rather extended distribution along 
the equator, whose normalized second moment is the minimum value 2/3. This result shows that 
|1, l) n can not be written as a tensor product of single-particle states. 

which are localized around the north pole (9 = tt), the equator (9 = vr/2) and the south pole 
(9 = 0) respectively. (See Fig. 0) Among the three states, \2, 0} n and |0,2) n are separable, 
but 

|1, l) n = -±=(|1, 0) n |0, l) n + |0, l) n |l, 0) n ) (47) 

is not. Coresponding to this fact, the Husimi distribution of |1, l) n is broader than the other 
two. 

C. moments 

The q-th moment of this distribution is denned as 

M$ = (2j + 1) j dfi(C) {n M (()} q , (48) 
where d/j,(£) is the Haar measure 




is invariant under the SU (2) transformation by definition of the Haar measure. 
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We can calculate the moments according to the general formula ()24|) and (|25|) . Here we 
consider the second moment, which can be calculated from the CG (Clebsch-Gordan) series 

An/2 ® D m/2 =D m @ An-1 © • • • A)- (50) 

If the state \<p) is represented as 

m2/ 

Iv) = c u\m/2,u) s , (51) 

f=— m/2 

the normalized second moment is 

2m + 1 



^ = J2 (52) 



m + 

H=—m 



= (m, fi\m/2, m/2; v, v')c v c v i. (53) 
Here, (m, /j l \m/2,m/2; u, v') is the CG coefficient, whose explicit form is 



(54) 



Let us examine some simple cases. If m = 1 (j = 1/2), 

<MA> = (|c 1/2 | 2 + |c_ 1/2 | 2 ) 2 = 1, (55) 

where P,- denotes the projection operator to the spin j representation. This result is natural 
because this is a one-body system and hence there is no many-body correlation. In the 
simplest non-trivial case m = 2 (j = 1), the normalized second moment is 

M (2) = (P2) = |ci| 4 + |c_i| 4 + 2|cic | 2 + 2|c c_i| 2 + ||cic_i + c 2 ,] 2 . (56) 

We can also write it as 

M(2) = ! _ (Pl) _ (Po) (57) 

according to (p?2*|) . (Pi) vanishes identically because the spin 1 part is anti-symmetric with 
respect to the exchange of the two components of the tensor product. The spin part is 

|0, 0), = ^=(\h l) s |l, -l) s - |1, 0) s |l, 0) s + |1, -l) s |l, 1),) . (58) 
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Hence 

M^ = l- l -\cl-2c lC ^\ 2 . (59) 

The maximum value of is 1, which is obtained when Cg = 2cic_i. This is the 
condition for the coherent states. The minimum value is 2/3, which is obtained if and only 
if 

3 {<?A<?-x) = 0, (60) 
3? (c^clj < 0, (61) 

N = |c_i|. (62) 
For example, |1, 1)„ satisfies this condition. 

IV. GENERAL TREATMENT OF BOSONIC SYSTEMS 

In this section we treat bosonic systems generally. 

A. Coherent states for bosonic systems 

Suppose there are m bosons in N single-particle states. The particle-hole excitation 
operators 

Xi = (na] (l<i,j<N) (63) 

satisfy the commutation relations of u(N): 

[Xi,Xf\=Xf8i-X*8{. (64) 

These operators generates the single-particle transformation: 

[Xi,a k } = 5{ ai . (65) 

The N single-particle states form the fundamental representation of U(N), and the 
bosonic many-body states form an irreducible representation which consists of symmet- 
ric combinations of the single-particle states. In the Young diagram, this representation is 
[l m ] (See Fig. EJ), and we denote it as B[N, m]. We will drop [N, m] when it is obvious. 
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qui 



[1**] = fl,l, ,1] 

B^[N,m] = B[N,qm] 

FIG. 3: Irreducible representations of U(N), B[N,m] = B^[N,m] and B^[N,m\. 

A coherent state is obtained by applying the U(N) transformation to the bosonic "vac- 
uum" |0) = |0, . . . ,0,m). As we have explained in the introduction, the bosonic coherent 
states are of the form 

10 = MOMQ) ••X0> (66) 

where \(f(C)) = U(()\<pi). Unitary transformations in the subspace spanned by |<^ 2 ), • • • , \<Pn) 
is irrelevant for the definition of the coherent state ()66jh and the phase transformation of 
\<Pi) is also irrelevant. Therefore the coherent states are specified by a point on the manifold 

U{N)/ (U{1) x U(N - 1)) ~ CP N -\ (67) 

An explicit form of the parametrization of the coherent states is 

|C)=Ar(C)exp^OX^ |0>, (68) 

where ( = (£ 2 , (3, ■ ■ ■ , Cat) is a coordinate system of CP^ -1 . 

The stationary phase condition of the coherent state path integral leads to an equation of 
motion on the phase space CP^™ 1 , which describes the dynamics of the condensed bosons. 

B. Moments of the Husimi distribution 

The Husimi distribution is defined as 

n M (0 = \(c\<p)\ 2 . (69) 

The normalized g-th moment thereof is 

^iS = l^pf M S = (dim ^ A) / Mc) i n \MY • (70) 
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Here, -Da is the bosonic Hilbert space B[N,m] = [l m ], and D q \ corresponds to [l qm ], which 
we denote as B^[N,m). (See Fig. EJ) Its dimension is 

dimgW[ jy>m] = ^ + ^- 1 )i , (71) 
L J (N-l)\(qm)\ v ' 



hence 

dimD qX _ r(iV + qm)T(m + 1) 



(72) 



dimD A T(N + m)T(qm + 1)' 
We can calculate the normalized moments by using the general formulae (}2"4*j) and (|25|). 
Let us consider the simplest case q = 2. The basis states of can be constructed as 

\B®,k) = Cn{l (2) if>) m ® |0) m . (73) 

Here, fc = . . . , fcjv) represents the occupation numbers, where the total number \k\ = 

2~2jkj is constrained to 2m. C is a normalization constant, |0) m = |0,...,0, m) is the 
m-particle vacuum, and 

= Xl®I + I®X{ (74) 
is a generator of the Lie algebra in B® 2 . By expanding ()73|h we obtain 

i^ 2 ),fc) = ^njE ^(^r^^^lio)™®^,,, (75) 



j=2 | n,=0 V 3 
ki k N f N 



\ / h \ I 

° e ■ • • e n u. ( x i) nj ® ( x i) fc ^ nj ® 

n 2 =0 n,iv=0 lj=2 ^ > 



(76) 



where 



N I n l 



i=2 

1 m! 



i — r- - / x \ n N fx \ n N—l / ±\ n 2 \ 1 / 

l[(Xip\0) m = (oj,) (aj^) ...(4) ai N+ - +n2A -^f-|0)o (77) 

,11(4) J |0>o (78) 



ml ntl . 



n!n). 79 

Til! 

Here, n = (n l5 . . . , n^r) represents the occupation numbers, where ri\ = m — (n^ + ■ ■ ■ + n 2 ). 
The factorial of the vector n means 

jv 

n\ = Y[ nj \. (80) 

3=1 
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By substituting (|79j) into (|75|). we obtain 



|£( 2 U> = C" V J— | n )®|fc-n>. 
^ V n (fe - n) 

|n|=m V v ; 

The normalization constant C" is determined by 

1 X ^ k\ 

~ 2s n\(k-n)\ 

\n\=m v ' 

(2m)! 
(m!) 2 



Hence 



| B W fc ) = JW y / fc! - | w )(8)lfc-n). 
The cases with larger q can be treated in the same way. The result is 

where 

A; \ J n , n ^i ! n i (if ni + n 2 H hn, = fe) 



j^na,... n g / I o (else) 

Suppose a state |<£>) is represented in the occupation number representation: 

\V) = Y^ c k\ k )- 

n 

From the general formulae (|24|) and (|25jh we obtain 

! fci 



[fe|= 



:?/(</ 



,(m!)9 ^ / fc V /2 
v y; |m|=-=|n 9 |=m v ' ' ' q 



V. FERMIONS 



In this section we study fermionic systems. 
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FIG. 4: Schematic picture of many-body Husimi distributions in a fermionic system. The localized 
wave packets (coherent states) represent independent-particle states (Slater determinants), and 
broad states represent correlated states. All Slater determinants are connected by the U(N) group, 
and treated equally. The TDHF equation defines a flow in this phase space. 

A. Coherent states for fermionic systems 

We consider m fermions in N single-particle states. The particle-hole excitation operators 



satisfy the same equations ()64j) and (jfj5)l as in the bosonic case. Hence they generate the 
single-particle transformation group U(N). Fermionic many-particle states form an irre- 
ducible representation which consists of anti-symmetric combinations of the single-particle 
states. In the Young diagram, this representation is denoted by [m] (See Fig. EJ), and we 
refer to it as F [N, m] . 

As we have explained in the introduction, the fermionic coherent states are Slater deter- 
minants, which can be written as 



where |y?j(C)) = U(()\<pi). Unitary transformations in the subspaces spanned by \ipx) ■ ■ ■ |<^ m ) 
and \ip m+ i) . . . \ (Pn) are irrelevant when we consider the coherent states. Therefore the 
coherent states are specified by a point on the complex Grassmann manifold 




(90) 



i ? ,0 = ^(WC)>M0>...|MC)» 



(91) 




(92) 
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[m] 



[m q ] = [m, m, . . . ,m] 



F[N, m] 



F^) [N, m] 



FIG. 5: Irreducible representations of U(N), F[N,m] = F^[N,m] and F^[N,m\. 
An explicit parametrization of the coherent states is given by [l| 



where |0) = \ 0, ... ,0,1, ... ,1) is the fermionic "vacuum", and ( = {Q} is a coordinate 
system of G m> ^(C). 

The "classical" equation of motion defined on the manifold is the TDHF (time-dependent 
Hartree-Fock) equation Q]. (See Fig. HJ) Therefore the manifold G mi Ar(C) is sometimes 
referred to as the TDHF manifold. The Husimi distribution on this manifold is defined as 




(93) 



N-m 



rn 



H W) {Q = \{F,t\ip)\ 2 . 



(94) 



B. Moments of the Husimi distribution 



Let us consider a m-particle state 




(95) 



n 



where 



n = (m, n 2 ,..., n N ) (m = or 1) 



(96) 



and 




(97) 
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We calculate the normalized g-th moment of the Husimi distribution 

#) = ^ M ft (98) 
1^ dim^A \f) y 1 

Here, D q \ is [m q ] in the Young diagram, which we refer to as F^[N,m] (Fig. EJ). In this 

case, 

f}„ (N-j-iy.(q + j)V 



hence 



dim£) gA " ' r(iV + g-j)r(j + 2) 

dimD A 11 T(iV+l-j)r(g + j + l)' 1 J 

According to the general formulae (f2"4*|) and (|25)1. 

%r = ^(F (<? U,r|(|n 1 )®---® |n,)) cn^'-cnq- (102) 
n 

Here an additional index r is introduced because the occupation number k does not neces- 
sarily specify a state in F™. The group-theoretical factors in ()102|) can be evaluated, for 
instance, by the eigenfunction method [ll| . However, it is not so easy as in the bosonic case 
to obtain an explicit formula for them. In the following, we calculate ()102)) in some simple 
examples. 



C. Examples 



Let us consider some simple fermionic systems. The cases with m — 1 are trivial, where 
is always equal to unity. Let us consider the cases with m = 2. It is obviously trivial 
for N = 2. For N = 3, there are three many-body states |011), 1 101) and |110), which form 
the representation 3 of U(3) (See Fig. |U1) The representation is 6, whose components 
are written as, for instance, 

1 1 n \ 

(103) 
(104) 
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F[3, 1] 

|010) 



Hi 



F[3,2] 



1 100) 



|110) 



H 



|001) 



|011) 



1 101) 



F( 2 )[3,2] |22 °W 2 



1121) 



|022) | 112) |202) 



1211) 



FIG. 6: Irreducible representations of U(3). The weight vector {H\,H<i) is related to occupation 
numbers by H\ = (rei — n<z)/2 and #2 = (^i + n 2 — 2n%) / (2y/?>) . Single-particle states form the 
fundamental representation F[3, 1] = 3, and two-particle (i.e., one-hole) states form F[3, 2] = 3. 
F( 2 )[3, 2] corresponds to 6. 



Here we arranged the elements of the tensor products vertically. For example, 



110 
101 



1 110> <g> |101). 



(105) 



Note that the sum of the two row vectors gives the occupation number of the product 

All the components of 6 are obtained by permutating the occupation numbers in (JI 
and (|l()4jl . Then we obtain the components of as 



B220 
-B211 



(ci 2 ) 2 , 
V^Cia C13, 



(106) 
(107) 



Here, Cy is the coefficient of the basis state where z-th and j-th states are occupied. For 
example, c\i is the coefficient of |011). The explicit formula for the second moment is 

m (2) = (M 2 + M 2 + |c 3 i| 2 ) 2 = i. (108) 
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^3 = ^ 

1 1210) 



v/6 12200) 
— H 2 



H 3 = -f 



Ho 



|2110) 



Mi 



|0202) | 1102) |2002) 
H x 



|0220) |1120) |2020) 



|0112) 



|1012) 



#3 = 



11201) 



H 2 

12101) 



- 0211 — 1111 — 0121 - 



|0022) 



Hi 



10121) 



11021) 



FIG. 7: An irreducible representation of f/(4), F( 2 )[4,2] = 20. The weight vector (Hi,H 2 ,H 3 ) 
is related to occupation numbers by H\ = (rti — n 2 )/2, H 2 = (rti + n2 — 2n3)/(2\/3) and #3 = 
(ni+n2+ri3 — 3n4)/(2\/6). This representation forms a regular octahedron in the three-dimensional 
weight space. 1 1111) is at the center thereof, and doubly degenerate. 

This result is natural because F[3, 2] is a one-hole system, and has no many-body corre- 
lation essentially. By the same reason, F[N, N — 1] is also trivial for arbitrary N. In this 
one-hole case, a state can be written as 



N 



(109) 



Here, \j) is a state with a hole in the j-th single-particle state. Since this system is essentially 
the same as the one boson system B[N, 1], we can use the general formulae for bosons 
Then, by using (|2*Hj) we obtain 



B 



k 



k\\ k 2 \ . . . fcjy! 



n ki fc 2 fcjv 
c l <-2 • • • C N 



(110) 



and 



in) 



The simplest non-trivial example is F[4, 2], whose dimension is six. To calculate the 
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second moment, we have to consider the CG series 



[2] ® [2] = [2,2] © [3,1] © [4] 
6 © 6 = 20 © 15 © 1. 

The components of = [2, 2] are written as, for instance, 

|2200) = 



;ii2) 



1100 
1100 



|2110) 



1 

7! 



1010 
1100 




(113) 
(114) 



Note that 1 1 1 1 1) is doubly degenerate. (See Fig. 0) We define two basis vectors as 
|llll,a) 

|1111,6> 





1100 \ 


0011 \ 


1010 \ 


0101 \ | 


j 












0011 / 


1100 / 


0101 / 


1010 / j 



(115) 



2^3 



-2 



J 


1100 \ + 


0011 \ 


1010 \ 


0101 




0011 / 


1100 / 


0101 / 


1010 



1001 \ 


0110 \ | 


)- 2 




0110 / 


1001 / j 



(116) 



Then 



-B2200 
-B2110 

Bmih 



(ci 2 ) 2 , 

V2c 12 C 13 , 

C12C34 + C13C24, 

C12C34 - C13C24 - 2ci 4 c 23 



(117) 
(118) 
(119) 
(120) 



The other components are obtained by permutating the occupation numbers. Then we 
obtain 



I 1 4 1 I 1 4 1 I 1 4 1 I 1 4 j I 14 1 I 4 
I C12 1 T I C13 1 -+- |Ci4| + I C23 1 + I c 24| + I C34 1 

+2 { IC12C13 | 2 + |Ci2Ci4| 2 + IC13C14 1 2 + IC12C23 | 2 + \C\2C24\ 2 + IC23C24I' 
+ |C13C 2 3| 2 + |ci 3 C 34 | 2 + |C 23 C 3 4| 2 + |ci 4 C 24 | 2 + |ci 4 C 34 | 2 + |C 24 C34| 2 } 
+ jj {IC12C34 + Ci3C 24 | 2 + |ci 2 C 34 - Ci4C 23 | 2 + |ci 3 C 2 4 + Ci4C 23 | 2 } 



(121) 
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A simpler expression of the second moment can be obtained from ()32j) : 



= 1-(P [3A] }-(P [A] ). 



(122) 



(-P[3,i]) vanishes identically because it is anti-symmetric with respect to the exchange of the 
components of the tensor product. Therefore all we have to calculate is (-P[4j). The basis 
state of the singlet representation [4] is 



l|4|> = 7E 



Hence 



1100 
0011 



0011 
1100 



1010 
0101 



0101 
1010 



1001 
0110 



0110 
1001 



(123) 



2 2 

1 — \P[4]) — 1 ~~ g IC12C34 — C13C24 + C14C23I 



(124) 



The second moment of the general two-fermion case F[N, 2] (N > 4) is obtained by a 
similar calculation. The result is 



M (2) = 5^|c n , i2 | 4 + 2 \ c Kh c Kh\ 2 + A (h,h,h,k) 



(125) 



Here, 



A(i,j,k,l) = - {|^-c fc , + c ifc c„l 2 + | Ci ,-c fc , - cc^l 2 + Ic,^ + cc^l 2 } , (126) 



and the sum is taken over 1 < %\ < < N, 1 < ji < ji < N, 1 < k < N and 1 < l\ < I2 < 
h < h < N. This can be rewritten as 



(127) 



h,hMM 

by using the constraint |cjii 2 | 2 = 1- 

VI. DISTINGUISHABLE PARTICLES 
A. General formalism 

Although all elementary particles are considered to be bosons or fermions, we often treat 
particles as distinguishable if the exchange effect is negligible. When we treat distinguish- 
able particles, we can change the single-particle basis of each particle separately. If there 
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are m particles and the z-th particle can take Ni states, the group of the single-particle 
transformation is 

SU(Nx) x SU(N 2 ) x • • • x SU(N m ). (128) 
We denote the Hilbert space of this system as D[Nx, N 2 , . . . , N m ], whose dimension is 

m 

m ' * * 

rii=i N{. It can also be written as [1] ® [1] <g> • • • ® [1], where z-th [1] is the fundamental rep- 
resentation of SU(Ni). A coherent states of this system is a separable (disentangled) state, 
which is considered to be a point on the "phase space" CP^ 1-1 x CP^ 2 " 1 x ■ • • x <CP Nm ^ 1 . 
The Husimi distribution is defined as a distribution function on this phase space. 



B. Moments 

The normalized g-th moment is defined as 

m = dimDqX M^. (129) 
M dimP A \f) V ; 

Here, D qX = [l q }<%>- ■ -<g> [l q ], which we denote as D^[Ni, . . . , N m }. We will drop [N h . . . , N m ] 

when it is obvious. 

Let us calculate the normalized moment by using the general formulae ()24|) and (|25|). If 
m — 1, this is the same as the bosonic case with m — 1. A state in [l q ] is represented by 
occupation numbers k = (hi, . . . , fcjv x ), which satisfies |fe| = g. Then 

fc|(|n 1 )®...®|n 9 )) = 4rf * V- ( 130 ) 
Here, |P^,fc) G [l 9 ] and |n;) G [1]. Since nj = 1, (HUJ) can be written simply as 

(DM,k\(\n 1 )®..-®\n q )) = < V q[ V J . (131) 

[ (else) 

The coefficients in (|25|) for multi-particle cases are easily obtained from (|131j) . A basis 
state in D^ q > [Ni, . . . , N m ] is specified by a set of occupation numbers {K,^}-, where 1 < z < m, 
1 < ji < N and J^-. kij. = q for Vz. Then we have 



Xh f^ 1 (if £,=i ";.;../ = for V(z, ./,) 



(^ (9) ,{^}|(|{^,1}>®---®|{W»= V V <1?2) 

[ (else) 

where | £><«>, {k^}) G D^^, ...,N m ) and IK,-,,,}) G P[iV\, . . . , N m ). 
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C. Qubits 



In this and the next subsections, we investigate systems of qubits rather in detail. Let 

m 

us consider a system of m-qubits, which is Dp., 2,"?. . , 2] in our notation. This is the spin 
(1/2, 1/2,..., 1/2) representation of SU(2) x SU(2) x ••• x SU(2), and D^[2, 2, . . . , 2] is 
the spin (q/2,q/2, . . . ,q/2) representation thereof. We denote a basis state of this system 
as |/Lt) = |/ii,/i 2 , . . . ,/i m ), where /ij takes 1/2 or —1/2. This notation is better than using 
1 and in that the symmetry between the two states is obvious. In the same way, a basis 
state in D^[2,2, . . . , 2] is denoted as |//i,// 2 , • • • , /-O, where runs from — q/2 to g/2 with 
the unity steps. This is related to the occupation number representation in the previous 
subsection by 

n it i = | + (133) 
n it2 = | - Hi- (134) 

According to the results in the previous subsection, the normalized g-th moment 

Khi^YK) ( 135 ) 

is obtained as 

^£ ) = Ei^i 2 ' ( 136 ) 

B v = {D^\u\{\^)®\^ 2 )®---®\^ q ))c^...c^ q (137) 

Til,..., Tig 



and 



(D® u\ (|^> | M2 > ® • • • |/* 9 » = { Utl { * 9+ J V2 ^ EJ=1 • (138) 

[ (else) 

D. Examples 

Let us start with a single qubit case D{2\. A qubit state is represented as 

\<p) = c + \+) + c-\-), (139) 
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where + and — are abbreviations of +1/2 and —1/2. The normalized q-th moment of the 
state is 

q/2 

m,S = £ i^i 2 > ( 14 °) 

where is obtained from (|2~%1) and (11381) as 



\ 1/2 

/J„ ( ^ J c + ql2+u cJ /2 - u . (141) 
Therefore 

M (?) = (\c + \ 2 + \c^\ 2 ) q = 1. (142) 

This case is trivial, as it should be. 

Next we consider 2-qubit case _D[2,2]. A state in this space is represented as 

\<p) = c ++ \ + +> + c+_| + -> + c_+| - +) + c __| - -). (143) 

D^[2,2] is a 9-dimensional representation, whose basis states are specified by v = (^1,^2) 
(— 1 < z/j < 1). The coefficients -Bi/ are easily obtained as 

Bu = (c ++ ) 2 , (144) 
B w = V2c ++ c + _, (145) 
B 00 = c ++ c__ + c+_c_+, (146) 



The other components which are not explicitly written here are obtained by exchanging the 
first and the second suffixes and changing their signs (+ <-» — , 1 <-> —1). 
The explicit form of the normalized second moment is 

M,S = |c ++ | 4 +|c + _| 4 + |c_ + | 4 + |c__| 4 

+2|c ++ c+_| 2 + 2|c ++ c h | 2 + 2|c Ch | 2 + 2|c c_ + | 2 

+ |c++c— + c + _c_+| 2 . (147) 

Let us rewrite this by using (J3*2~j) as 

M$ = 1 - (P (li0) > - (P ( o,i)) - (P ( o,o)), (148) 
26 



where ( Ji, J 2) represents the irreducible representation with the total spins J\ and J 2- The 
spin 1 representation, which has three components 



|1> 



+ \ , x 1 

|0) = 7I 




1} 



(149) 



are symmetric with respect to the exchange between the components of the tensor product. 
The spin representation, whose basis state is 



10} 



1 



+ 



(150) 



is anti-symmetric with respect to the exchange. Note that we arranged the components of the 
tensor product vertically so that we can distinguish it from the tensor product to represent 
the physical composite systems. (-P(i,o)) an d (-P(o,i)) van ish because the representations (1, 0) 
and (0, 1) are anti-symmetric as a whole, and only (P(o,o)) survives in ()148j) . The basis state 
of the singlet representation (0, 0) is 



|00) = 7=2 




(151) 
(152) 



Then we obtain 



C 2 

(P(o,o)) = \c++c— - c+_c_+| 2 = — , 



(153) 



where C is the concurrence 



121. Hence 



m = 1 - — . 



The next example is the second moment of the 3-qubit case. 



(154) 



-Bin 
-Bno 
-B100 
-Bqoo 



(c +++ ) 2 , 



y/2, 



C_|_- L. |_C_I 



•+-) 



C +++ C H + c ++- c +- + ) 

1 , 



V2 



Ci 1 1 c_ 



+ c 4 



+ c_++c+_). 



(155) 
(156) 
(157) 
(158) 
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The other components are obtained by permutating the suffixes and changing their signs. 
Then we obtain 



M, 



(2) 
\<P) 



|C+H 



+ C. 



+ |c ++ _| 4 + |c H 



| 4 + |c_ ++ | 4 +|c__ + | 4 +|c. 



+ |c+- 



'2 c I i I cj 



2|c +++ c+_ + | 2 + 2|c +++ c_ ++ | 2 + 2|c___c_ 



+2|c c_ + _| 2 + 2|c c + | 2 + 2|c + _ + c+__| 2 + 2|c_ ++ c + | 2 

+2|c ++ _c_ + _| 2 + 2|c_ ++ c_+_| 2 + 2|c++_c + __| 2 + 2|c + _ + c +| 2 



+ I c +++ C H + c ++- c +-+\ + l c +++ c 



c+_+c__ 



I 2 + l c +++ c -+- + c -++ c ++-l 2 



+|c c_ ++ + c — + c_ + _| 2 + |c c ++ _ + c_ + _c + — | 2 + |c C + _ + + C + c — + | 2 



C++_C_ 



C+- + C- 



c -++ c +- 



(159) 



By using (|32l), can also be written as 



= i - (p (0A1) > - (P (0 ,i,o)> - (P(im)- 



(160) 



Here we dropped the anti-symmetric representations (1, 1, 0), (1, 0, 1), (0, 1, 1) and (0, 0, 0). 
(0, 0, 1) is a 3-dimensional representation, whose basis states are 



|001) 



|000) 



+ + 



+ + + 



+ - + 




(161) 



2V2 



+ 



+ 



+ 



+ 



+ + 



+ + + 




100-1) 



+ 



(162) 



(163) 



Hence 



1 

t:I c +++ c — h — c +-+ c -- 



+ 2l c +++ c - 



| 2 + |c ++ _c___ - c + __c_ + _| 2 

|2 



c+-+c- 



C_| C- 



(164) 
(165) 



By explicit calculation this is shown to be equal to ^Tr (pabPab)- Here, A,B and C are the 
three qubits. pab is the density matrix of the pair A and B which is obtained by tracing 
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out the information of C. pab is the "spin-flipped" density matrix 



1131 



PAB = (<Ty ® Vy)p*AB(. a V ® <fy)j 



(166) 



where the asterisk denotes complex conjugation in the standard basis and 








i 



— % 



) 



(167) 



According to [1 



Tr (pabPab) — Cab + 7^ t abc, 



(168) 



where Cab is the concurrence of the density matrix pab and tabc is the 3-t angle whose 
explicit form is given in [ljj]. In the same way we can calculate (P(o,i,o)) an d (P(o,o,i)), and 
obtain the result 



VII. CONCLUDING REMARKS 

In this paper, we have defined generalized Husimi distributions for many-body systems, 
and showed that their moments can be practically useful measures of complexity of many- 
body quantum states. We have derived an algebraic formula to evaluate the moments 
with integer indices. For bosons and distinguishable particles we have derived the group 
theoretical factors appearing in the formula explicitly. We have also examined some special 
cases of fermions, though we have not yet obtained a general formula for the group-theoretical 
factors in fermionic cases. 

The moments take the maximum values if and only if the state is an independent-particle 
state. Therefore the moments are considered to be measures of entanglement. Although a 
lot of measures of entanglement have been proposed so far, our method has an advantage of 
being able to produce infinitely many measures systematically. Of course, not all of them are 
independent. Nevertheless, it might be possible that the moments and the related invariants 
(See f)32jl ) form a complete set of algebraic invariants to classify the pure state entanglement. 

The relation between derealization of the Husimi distribution and chaoticity of the clas- 
sical mechanics has been shown in one-body systems We can expect a similar 




1 



3 



(169) 
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correspondence also in many-body systems, which means a relation between many-body 
correlation and chaoticity of the dynamics of coherent states (mean fields). For example, 
chaotic behaviors in fermionic systems like nuclei have been studied from the viewpoints of 

n n 

the shell model [16|] and the mean field dynamics [10J. However, the relation between the 
two standpoints has been unclear. We hope that the Husimi distribution will be a bridge 
between the two fields. 

In this paper, we have concerned only with the systems with definite particle numbers, 
and considered classifications based on the single-particle unitary transformations. When 
we treat a systems with variable particle number, we can introduce a generalized one-body 
transformation, which is known as the Bogoliubov transformation. Then, for example, BCS- 
type wavefunctions are considered to be coherent states in fermionic systems. The treatment 
of such systems will be reported elsewhere. 
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